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PREFACE 



This work began as a collaborative attempt with David Halliday to revise and 
update the second edition of his classic text Introductory Nuclear Physics (New 
York: Wiley, 1955). As the project evolved, it became clear that, owing to other 
commitments, Professor Halliday would be able to devote only limited time to 
the project and he therefore volunteered to remove himself from active participa- 
tion, a proposal to which I reluctantly and regretfully agreed. He was kind 
enough to sign over to me the rights to use the material from the previous edition. 

I first encountered Halliday’s text as an undergraduate physics major, and it 
was perhaps my first real introduction to nuclear physics. I recall being impressed 
by its clarity and its readability, and in preparing this new version, I have tried to 
preserve these elements, which are among the strengths of the previous work. 

Audience This text is written primarily for an undergraduate audience, but 
could be used in introductory graduate surveys of nuclear physics as well. It can 
be used specifically for physics majors as part of a survey of modern physics, but 
could (with an appropriate selection of material) serve as an introductory course 
for other areas of nuclear science and technology, including nuclear chemistry, 
nuclear engineering, radiation biology, and nuclear medicine. 

Background It is expected that students have a previous background in quan- 
tum physics, either at the introductory level [such as the author’s text Modern 
Physics (New York: Wiley, 1983)] or at a more advanced, but still undergraduate 
level. (A brief summary of the needed quantum background is given in Chapter 
2.) The text is therefore designed in a “ two-track” mode, so that the material that 
requires the advanced work in quantum mechanics, for instance, transition 
probabilities or matrix elements, can be separated from the rest of the text by 
skipping those sections that require such a background. This can be done without 
interrupting the logical flow of the discussion. 

Mathematical background at the level of differential equations should be 
sufficient for most applications. 

Emphasis There are two features that distinguish the present book. The first is 
the emphasis on breadth. The presentation of a broad selection of material 
permits the instructor to tailor a curriculum to meet the needs of any particular 
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student audience. The complete text is somewhat short for a full-year course, but 
too long for a course of quarter or semester length. The instructor is therefore 
able to select material that will provide students with the broadest possible 
introduction to the field of nuclear physics, consistent with the time available for 
the course. 

The second feature is the unabashedly experimental and phenomenological 
emphasis and orientation of the presentation. The discussions of decay and 
reaction phenomena are accompanied with examples of experimental studies 
from the literature. These examples have been carefully selected following 
searches for papers that present data in the clearest possible manner and that 
relate most directly to the matter under discussion. These original experiments 
are discussed, often with accompanying diagrams of apparatus, and results with 
uncertainties are given, all in the attempt to convince students that progress in 
nuclear physics sprang not exclusively from the forehead of Fermi, but instead 
has been painstakingly won in the laboratory. At the same time, the rationale and 
motivation for the experiments are discussed, and their contributions to the 
theory are emphasized. 

Organization The book is divided into four units: Basic Nuclear Structure, 
Nuclear Decay and Radioactivity, Nuclear Reactions, and Extensions and Appli- 
cations. The first unit presents background material on nuclear sizes and shapes, 
discusses the two-nucleon problem, and presents an introduction to nuclear 
models. These latter two topics can be skipped without loss of continuity in an 
abbreviated course. The second unit on decay and radioactivity presents the 
traditional topics, with new material included to bring nuclear decay nearly into 
the current era (the recently discovered “heavy” decay modes, such as 14 C, 
double fi decay, /^-delayed nucleon emission, Mossbauer effect, and so on). The 
third unit surveys nuclear reactions, including fission and fusion and their 
applications. The final unit deals with topics that fall only loosely under the 
nuclear physics classification, including hyperfine interactions, particle physics, 
nuclear astrophysics, and general applications including nuclear medicine. The 
emphasis here is on the overlap with other physics and nonphysics specialties, 
including atomic physics, high-energy physics, cosmology, chemistry, and medi- 
cine. Much of this material, particularly in Chapters 18 and 19, represents 
accomplishments of the last couple of years and therefore, as in all such volatile 
areas, may be outdated before the book is published. Even if this should occur, 
however, the instructor is presented with a golden opportunity to make important 
points about progress in science. Chapter 20 features applications involving 
similarly recent developments, such as PET scans. The material in this last unit 
builds to a considerable degree on the previous material; it would be very unwise, 
for example, to attempt the material on meson physics or particle physics without 
a firm grounding in nuclear reactions. 

Sequence Chapters or sections that can be omitted without loss of continuity in 
an abbreviated reading are indicated with asterisks (*) in the table of contents. 
An introductory short course in nuclear physics could be based on Chapters 1, 2, 
3, 6, 8, 9, 10, and 11, which cover the fundamental aspects of nuclear decay and 
reactions, but little of nuclear structure. Fission and fusion can be added from 
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Chapters 13 and 14. Detectors and accelerators can be included with material 
selected from Chapters 7 and 15. 

The last unit (Chapters 16 to 20) deals with applications and does not 
necessarily follow Chapter 15 in sequence. In fact, most of this material could be 
incorporated at any time after Chapter 11 (Nuclear Reactions). Chapter 16, 
covering spins and moments, could even be moved into the first unit after 
Chapter 3. Chapter 19 (Nuclear Astrophysics) requires background material on 
fission and fusion from Chapters 13 and 14. 

Most of the text can be understood with only a minimal background in 
quantum mechanics. Chapters or sections that require a greater background (but 
still at the undergraduate level) are indicated in the table of contents with a 
dagger (f). 

Many undergraduates, in my experience, struggle with even the most basic 
aspects of the quantum theory of angular momentum, and more abstract con- 
cepts, such as isospin, can present them with serious difficulties. For this reason, 
the introduction of isospin is delayed until it is absolutely necessary in Chapter 
11 (Nuclear Reactions) where references to its application to beta and gamma 
decays are given to show its importance to those cases as well. No attempt is 
made to use isospin coupling theory to calculate amplitudes or cross sections. In 
an abbreviated coverage, it is therefore possible to omit completely any discus- 
sion of isospin, but it absolutely must be included before attempting Chapters 17 
and 18 on meson and particle physics. 

Notation Standard notation has been adopted, which unfortunately overworks 
the symbol T to represent kinetic energy, temperature, and isospin. The particle 
physicist’s choice of / for isospin and J for nuclear spin leaves no obvious 
alternative for the total electronic angular momentum. Therefore, I has been 
reserved for the total nuclear angular momentum, J for the total electronic 
angular momentum, and T for the isospin. To be consistent, the same scheme is 
extended into the particle physics regime in Chapters 17 and 18, even though it 
may be contrary to the generally accepted notation in particle physics. The 
lowercase j refers to the total angular momentum of a single nucleon or atomic 
electron. 

References No attempt has been made to produce an historically accurate set of 
references to original work. This omission is done partly out of my insecurity 
about assuming the role of historian of science and partly out of the conviction 
that references tend to clutter, rather than illuminate, textbooks that are aimed 
largely at undergraduates. Historical discussions have been kept to a minimum, 
although major insights are identified with their sources. The history of nuclear 
physics, which so closely accompanies the revolutions wrought in twentieth-cen- 
tury physics by relativity and quantum theory, is a fascinating study in itself, and 
I encourage serious students to pursue it. In stark contrast to modern works, the 
classic papers are surprisingly readable. Many references to these early papers 
can be found in Halliday’s book or in the collection by Robert T. Beyer, 
Foundations of Nuclear Physics (New York: Dover, 1949), which contains reprints 
of 13 pivotal papers and a classified bibliography of essentially every nuclear 
physics publication up to 1947. 
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Each chapter in this textbook is followed with a list of references for further 
reading, where more detailed or extensive treatments can be found. Included in 
the lists are review papers as well as popular-level books and articles. 

Several of the end-of-chapter problems require use of systematic tabulations of 
nuclear properties, for which the student should have ready access to the current 
edition of the Table of Isotopes or to a complete collection of the Nuclear Data 
Sheets . 
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BASIC 
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STRUCTURE 




I ll 

BASIC CONCEPTS 



Whether we date the origin of nuclear physics from Becquerel’s discovery of 
radioactivity in 1896 or Rutherford’s hypothesis of the existence of the nucleus in 
1911, it is clear that experimental and theoretical studies in nuclear physics have 
played a prominent role in the development of twentieth century physics. As a 
result of these studies, a chronology of which is given on the inside of the front 
cover of this book, we have today a reasonably good understanding of the 
properties of nuclei and of the structure that is responsible for those properties. 
Furthermore, techniques of nuclear physics have important applications in other 
areas, including atomic and solid-state physics. Laboratory experiments in nuclear 
physics have been applied to the understanding of an incredible variety of 
problems, from the interactions of quarks (the most fundamental particles of 
which matter is composed), to the processes that occurred during the early 
evolution of the universe just after the Big Bang. Today physicians use techniques 
learned from nuclear physics experiments to perform diagnosis and therapy in 
areas deep inside the body without recourse to surgery; but other techniques 
learned from nuclear physics experiments are used to build fearsome weapons of 
mass destruction, whose proliferation is a constant threat to our future. No other 
field of science comes readily to mind in which theory encompasses so broad a 
spectrum, from the most microscopic to the cosmic, nor is there another field in 
which direct applications of basic research contain the potential for the ultimate 
limits of good and evil. 

Nuclear physics lacks a coherent theoretical formulation that would permit us 
to analyze and interpret all phenomena in a fundamental way; atomic physics 
has such a formulation in quantum electrodynamics, which permits calculations 
of some observable quantities to more than six significant figures. As a result, we 
must discuss nuclear physics in a phenomenological way, using a different 
formulation to describe each different type of phenomenon, such as a decay, /? 
decay, direct reactions, or fission. Within each type, our ability to interpret 
experimental results and predict new results is relatively complete, yet the 
methods and formulation that apply to one phenomenon often are not applicable 
to another. In place of a single unifying theory there are islands of coherent 
knowledge in a sea of seemingly uncorrelated observations. Some of the most 
fundamental problems of nuclear physics, such as the exact nature of the forces 
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that hold the nucleus together, are yet unsolved. In recent years, much progress 
has been made toward understanding the basic force between the quarks that are 
the ultimate constituents of matter, and indeed attempts have been made at 
applying this knowledge to nuclei, but these efforts have thus far not contributed 
to the clarification of nuclear properties. 

We therefore adopt in this text the phenomenological approach, discussing 
each type of measurement, the theoretical formulation used in its analysis, and 
the insight into nuclear structure gained from its interpretation. We begin with a 
summary of the basic aspects of nuclear theory, and then turn to the experiments 
that contribute to our knowledge of structure, first radioactive decay and then 
nuclear reactions. Finally, we discuss special topics that contribute to micro- 
scopic nuclear structure, the relationship of nuclear physics to other disciplines, 
and applications to other areas of research and technology. 



1 . 1 HISTORY AND OVERVIEW 

The search for the fundamental nature of matter had its beginnings in the 
speculations of the early Greek philosophers; in particular, Democritus in the 
fourth century b.c. believed that each kind of material could be subdivided into 
smaller and smaller bits until one reached the very limit beyond which no further 
division was possible. This atom of material, invisible to the naked eye, was to 
Democritus the basic constituent particle of matter. For the next 2400 years, this 
idea remained only a speculation, until investigators in the early nineteenth 
century applied the methods of experimental science to this problem and from 
their studies obtained the evidence needed to raise the idea of atomism to the 
level of a full-fledged scientific theory. Today, with our tendency toward the 
specialization and compartmentalization of science, we would probably classify 
these early scientists (Dalton, Avogadro, Faraday) as chemists. Once the chemists 
had elucidated the kinds of atoms, the rules governing their combinations in 
matter, and their systematic classification (Mendeleev’s periodic table), it was 
only natural that the next step would be a study of the fundamental properties of 
individual atoms of the various elements, an activity that we would today classify 
as atomic physics. These studies led to the discovery in 1896 by Becquerel of the 
radioactivity of certain species of atoms and to the further identification of 
radioactive substances by the Curies in 1898. Rutherford next took up the study 
of these radiations and their properties; once he had achieved an understanding 
of the nature of the radiations, he turned them around and used them as probes 
of the atoms themselves. In the process he proposed in 1911 the existence of the 
atomic nucleus, the confirmation of which (through the painstaking experiments 
of Geiger and Marsden) provided a new branch of science, nuclear physics, 
dedicated to studying matter at its most fundamental level. Investigations into 
the properties of the nucleus have continued from Rutherford’s time to the 
present. In the 1940s and 1950s, it was discovered that there was yet another level 
of structure even more elementary and fundamental than the nucleus. Studies of 
the particles that contribute to the structure at this level are today carried out in 
the realm of elementary particle (or high energy) physics. 

Thus nuclear physics can be regarded as the descendent of chemistry and 
atomic physics and in turn the progenitor of particle physics. Although nuclear 
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physics no longer occupies center stage in the search for the ultimate components 
of matter, experiments with nuclei continue to contribute to the understanding of 
basic interactions. Investigation of nuclear properties and the laws governing the 
structure of nuclei is an active and productive area of physical research in its own 
right, and practical applications, such as smoke detectors, cardiac pacemakers, 
and medical imaging devices, have become common. Thus nuclear physics has in 
reality three aspects: probing the fundamental particles and their interactions, 
classifying and interpreting the properties of nuclei, and providing technological 
advances that benefit society. 



1 .2 SOME INTRODUCTORY TERMINOLOGY 

A nuclear species is characterized by the total amount of positive charge in the 
nucleus and by its total number of mass units. The net nuclear charge is equal to 
+ Ze, where Z is the atomic number and e is the magnitude of the electronic 
charge. The fundamental positively charged particle in the nucleus is the proton , 
which is the nucleus of the simplest atom, hydrogen. A nucleus of atomic number 
Z therefore contains Z protons, and an electrically neutral atom therefore must 
contain Z negatively charged electrons. Since the mass of the electrons is 
negligible compared with the proton mass (m p - 2000 m e ), the electron can often 
be ignored in discussions of the mass of an atom. The mass number of a nuclear 
species, indicated by the symbol A, is the integer nearest to the ratio between the 
nuclear mass and the fundamental mass unit, defined so that the proton has a 
mass of nearly one unit. (We will discuss mass units in more detail in Chapter 3.) 
For nearly all nuclei, A is greater than Z, in most cases by a factor of two or 
more. Thus there must be other massive components in the nucleus. Before 1932, 
it was believed that the nucleus contained A protons, in order to provide the 
proper mass, along with A - Z nuclear electrons to give a net positive charge of 
Ze . However, the presence of electrons within the nucleus is unsatisfactory for 
several reasons: 

1. The nuclear electrons would need to be bound to the protons by a very 
strong force, stronger even than the Coulomb force. Yet no evidence for this 
strong force exists between protons and atomic electrons. 

2. If we were to confine electrons in a region of space as small as a nucleus 
(Ax - 10 -14 m), the uncertainty principle would require that these electrons 
have a momentum distribution with a range A/? - h/tsx = 20 MeV/c. 
Electrons that are emitted from the nucleus in radioactive /? decay have 
energies generally less than 1 MeV; never do we see decay electrons with 
20 MeV energies. Thus the existence of 20 MeV electrons in the nucleus is 
not confirmed by observation. 

3. The total intrinsic angular momentum (spin) of nuclei for which A - Z is 
odd would disagree with observed values if A protons and A — Z electrons 
were present in the nucleus. Consider the nucleus of deuterium {A = 2, 
Z = 1), which according to the proton-electron hypothesis would contain 2 
protons and 1 electron. The proton and electron each have intrinsic angular 
momentum (spin) of and the quantum mechanical rules for adding spins 
of particles would require that these three spins of \ combine to a total of 
either § or Yet the observed spin of the deuterium nucleus is 1. 
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4. Nuclei containing unpaired electrons would be expected to have magnetic 
dipole moments far greater than those observed. If a single electron were 
present in a deuterium nucleus, for example, we would expect the nucleus to 
have a magnetic dipole moment about the same size as that of an electron, 
but the observed magnetic moment of the deuterium nucleus is about of 
the electron’s magnetic moment. 

Of course it is possible to invent all sorts of ad hoc reasons for the above 
arguments to be wrong, but the necessity for doing so was eliminated in 1932 
when the neutron was discovered by Chadwick. The neutron is electrically neutral 
and has a mass about equal to the proton mass (actually about 0.1% larger). Thus 
a nucleus with Z protons and A - Z neutrons has the proper total mass and 
charge, without the need to introduce nuclear electrons. When we wish to 
indicate a specific nuclear species, or nuclide , we generally use the form 
where X is the chemical symbol and N is the neutron number , A - Z. The 
symbols for some nuclides are }H 0 , 26 Fe 30 . The chemical symbol and the 

atomic number Z are redundant — every H nucleus has Z = 1, every U nucleus 
has Z = 92, and so on. It is therefore not necessary to write Z. It is also not 
necessary to write N, since we can always find it from A - Z. Thus 238 U is a 
perfectly valid way to indicate that particular nuclide; a glance at the periodic 
table tells us that U has Z = 92, and therefore 238 U has 238 — 92 = 146 
neutrons. You may find the symbols for nuclides written sometimes with Z and 
N, and sometimes without them. When we are trying to balance Z and TV in a 
decay or reaction process, it is convenient to have them written down; at other 
times it is cumbersome and unnecessary to write them. 

Neutrons and protons are the two members of the family of nucleons. When we 
wish simply to discuss nuclear particles without reference to whether they are 
protons or neutrons, we use the term nucleons. Thus a nucleus of mass number A 
contains A nucleons. 

When we analyze samples of many naturally occurring elements, we find that 
nuclides with a given atomic number can have several different mass numbers; 
that is, a nuclide with Z protons can have a variety of different neutron numbers. 
Nuclides with the same proton number but different neutron numbers are called 
isotopes ; for example, the element chlorine has two isotopes that are stable 
against radioactive decay, 35 Cl and 37 C1. It also has many other unstable isotopes 
that are artificially produced in nuclear reactions; these are the radioactive 
isotopes (or radioisotopes) of Cl. 

It is often convenient to refer to a sequence of nuclides with the same N but 
different Z; these are called isotones. The stable isotones with N = 1 are 2 H and 
3 He. Nuclides with the same mass number A are known as isobars ; thus stable 
3 He and radioactive 3 H are isobars. 



1 .3 NUCLEAR PROPERTIES 

Once we have identified a nuclide, we can then set about to measure its 
properties, among which (to be discussed later in this text) are mass, radius, 
relative abundance (for stable nuclides), decay modes and half-lives (for radioac- 
tive nuclides), reaction modes and cross sections, spin, magnetic dipole and 
electric quadrupole moments, and excited states. Thus far we have identified 
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Neutron number N 

Figure 1.1 Stable nuclei are shown in dark shading and known radioactive 
nuclei are in light shading. 



nuclides with 108 different atomic numbers (0 to 107); counting all the different 
isotope's, the total number of nuclides is well over 1000, and the number of 
carefully studied new nuclides is growing rapidly owing to new accelerators 
dedicated to studying the isotopes far from their stable isobars. Figure 1.1 shows 
a representation of the stable and known radioactive nuclides. 

As one might expect, cataloging all of the measured properties of these many 
nuclides is a formidable task. An equally formidable task is the retrieval of that 
information: if we require the best current experimental value of the decay modes 
of an isotope or the spin and magnetic moment of another, where do we look? 

Nuclear physicists generally publish the results of their investigations in 
journals that are read by other nuclear physicists; in this way, researchers from 
distant laboratories are aware of one another’s activities and can exchange ideas. 
Some of the more common journals in which to find such communications are 
Physical Review , Section C (abbreviated Phys. Rev. C), Physical Review Letters 
( Phys . Rev. Lett.), Physics Letters , Section B (Phys. Lett. B ), Nuclear Physics , 
Section A ( Nucl . Phys. A ), Zeitschrift fur Physik , Section A (Z. Phys. A ), and 
Journal of Physics , Section G ( J . Phys. G). These journals are generally published 
monthly, and by reading them (or by scanning the table of contents), we can find 
out about the results of different researchers. Many college and university 
libraries subscribe to these journals, and the study of nuclear physics is often 
aided by browsing through a selection of current research papers. 

Unfortunately, browsing through current journals usually does not help us to 
locate the specific nuclear physics information we are seeking, unless we happen 
to stumble across an article on that topic. For this reason, there are many sources 
of compiled nuclear physics information that summarize nuclear properties and 
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give references to the literature where the original publication may be consulted. 
A one-volume summary of the properties of all known nuclides is the Table of 
Isotopes , edited by M. Lederer and V. Shirley (New York: Wiley, 1978). A copy 
of this indispensible work is owned by every nuclear physicist. A more current 
updating of nuclear data can be found in the Nuclear Data Sheets , which not 
only publish regular updated collections of information for each set of isobars, 
but also give an annual summary of all published papers in nuclear physics, 
classified by nuclide. This information is published in journal form and is also 
carried by many libraries. It is therefore a relatively easy process to check the 
recently published work concerning a certain nuclide. 

Two other review works are the Atomic Data and Nuclear Data Tables , which 
regularly produces compilations of nuclear properties (for example, j8 or y 
transition rates or fission energies), and the Annual Review of Nuclear and Particle 
Science (formerly called the Annual Review of Nuclear Science ), which each year 
publishes a collection of review papers on current topics in nuclear and particle 
physics. 



1 .4 UNITS AND DIMENSIONS 

In nuclear physics we encounter lengths of the order of 10" 15 m, which is one 
femtometer (fm). This unit is colloquially known as one fermi, in honor of the 
pioneer Italian-American nuclear physicist, Enrico Fermi. Nuclear sizes range 
from about 1 fm for a single nucleon to about 7 fm for the heaviest nuclei. 

The time scale of nuclear phenomena has an enormous range. Some nuclei, 
such as 5 He or 8 Be, break apart in times of the order of 10" 20 s. Many nuclear 
reactions take place on this time scale, which is roughly the length of time that 
the reacting nuclei are within range of each other’s nuclear force. Electromagnetic 
(y) decays of nuclei occur generally within lifetimes of the order of 10" 9 s 
(nanosecond, ns) to 10" 12 s (picosecond, ps), but many decays occur with much 
shorter or longer lifetimes, a and /3 decays occur with even longer lifetimes, often 
minutes or hours, but sometimes thousands or even millions of years. 

Nuclear energies are conveniently measured in millions of electron-volts (MeV), 
where 1 eV = 1.602 X 10" 19 J is the energy gained by a single unit of electronic 
charge when accelerated through a potential difference of one volt. Typical ft and 
y decay energies are in the range of 1 MeV, and low-energy nuclear reactions take 
place with kinetic energies of order 10 MeV. Such energies are far smaller than 
the nuclear rest energies, and so we are justified in using nonrelativistic formulas 
for energy and momentum of the nucleons, but /3-decay electrons must be treated 
relativistically. 

Nuclear masses are measured in terms of the unified atomic mass unit , u, 
defined such that the mass of an atom of 12 C is exactly 12 u. Thus the nucleons 
have masses of approximately 1 u. In analyzing nuclear decays and reactions, we 
generally work with mass energies rather than with the masses themselves. The 
conversion factor is 1 u = 931.502 MeV, so the nucleons have mass energies of 
approximately 1000 MeV. The conversion of mass to energy is of course done 
using the fundamental result from special relativity, E = me 2 ; thus we are free to 
work either with masses or energies at our convenience, and in these units 
c 2 = 931.502 MeV/u. 
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REFERENCES FOR ADDITIONAL READING 

The following comprehensive nuclear physics texts provide explanations or 
formulations alternative to those of this book. Those at the introductory level are 
at about the same level as the present text; higher-level texts often form the basis 
for more advanced graduate courses in nuclear physics. No attempt has been 
made to produce a complete list of reference works; rather, these are the ones the 
author has found most useful in preparing this book. 

These “classic” texts now mostly outdated but still containing much useful 
material are interesting for gaining historical perspective: R. D. Evans, The 
Atomic Nucleus (New York: McGraw-Hill, 1955) (For 20 years, since his 
graduate-student days, the most frequently used book on the author’s shelves. Its 
binding has all but deteriorated, but its completeness and clarity remain.); David 
Halliday, Introductory Nuclear Physics (New York: Wiley, 1955); I. Kaplan, 
Nuclear Physics (Reading, MA: Addison-Wesley, 1955). 

Introductory texts complementary to this text are: W. E. Burcham, Nuclear 
Physics : An Introduction (London: Longman, 1973); B. L. Cohen, Concepts of 
Nuclear Physics (New York: McGraw-Hill, 1971); Harald A. Enge, Introduction 
to Nuclear Physics (Reading, MA: Addison-Wesley, 1966); Robert A. Howard, 
Nuclear Physics (Belmont, CA: Wadsworth, 1963); Walter E. Meyerhof, Ele- 
ments of Nuclear Physics (New York: McGraw-Hill, 1967); Haro Von Buttlar, 
Nuclear Physics : An Introduction (New York: Academic Press, 1968). 

Intermediate texts, covering much the same material as the present one but 
distinguished primarily by a more rigorous use of quantum mechanics, are: M. G. 
Bowler, Nuclear Physics (Oxford: Pergamon, 1973); Emilio Segre, Nuclei and 
Particles (Reading, MA: W. A. Benjamin, 1977). 

Advanced texts, primarily for graduate courses, but still containing much 
material of a more basic nature, are: Hans Frauenfelder and Ernest M. Henley, 
Subatomic Physics (Englewood Cliffs, NJ: Prentice-Hall, 1974); M. A. Preston, 
Physics of the Nucleus (Reading, MA: Addison-Wesley, 1962). 

Advanced works, more monographs than texts in nature, are: John M. Blatt 
and Victor F. Weisskopf, Theoretical Nuclear Physics (New York: Wiley, 1952); 
A. Bohr and B. R. Mottelson, Nuclear Structure (New York: W. A. Benjamin, 
1969); A. deShalit and H. Feshbach, Theoretical Nuclear Physics (New York: 
Wiley, 1974). 
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ELEMENTS OF 
QUANTUM MECHANICS 



Nucleons in a nucleus do not behave like classical particles, colliding like billiard 
balls. Instead, the wave behavior of the nucleons determines the properties of the 
nucleus, and to analyze this behavior requires that we use the mathematical 
techniques of quantum mechanics. 

From a variety of scattering experiments, we know that the nucleons in a 
nucleus are in motion with kinetic energies of the order of 10 MeV. This energy is 
small compared with the nucleon rest energy (about 1000 MeV), and so we can 
with confidence use nonrelativistic quantum mechanics. 

To give a complete introduction to quantum mechanics would require a text 
larger than the present one. In this chapter, we summarize some of the important 
concepts that we will need later in this book. We assume a previous introduction 
to the concepts of modern physics and a familiarity with some of the early 
experiments that could not be understood using classical physics; these experi- 
ments include thermal (blackbody) radiation, Compton scattering, and the photo- 
electric effect. At the end of this chapter is a list of several introductory modern 
physics texts for review. Included in the list are more advanced quantum physics 
texts, which contain more complete and rigorous discussions of the topics 
summarized in this chapter. 



2.1 QUANTUM BEHAVIOR 

Quantum mechanics is a mathematical formulation that enables us to calculate 
the wave behavior of material particles. It is not at all a priori evident that such 
behavior should occur, but the suggestion follows by analogy with the quantum 
behavior of light. Before 1900, light was generally believed to be a wave 
phenomenon, but the work of Planck in 1900 (analyzing blackbody radiation) 
and Einstein in 1905 (analyzing the photoelectric effect) showed that it was also 
necessary to consider light as if its energy were delivered not smoothly and 
continuously as a wave but instead in concentrated bundles or “quanta,” in effect 
“particles of light.” 

The analogy between matter and light was made in 1924 by de Broglie, 
drawing on the previous work of Einstein and Compton. If light, which we 
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generally regard as a wave phenomenon, also has particle aspects, then (so de 
Broglie argued) might not matter, which we generally regard as composed of 
particles, also have a wave aspect? Again proceeding by analogy with light, de 
Broglie postulated that associated with a “particle” moving with momentum p is 
a “wave” of wavelength X = h/p where h is Planck’s constant. The wavelength 
defined in this way is generally called the de Broglie wavelength. Experimental 
confirmation of de Broglie’s hypothesis soon followed in 1927 through the 
experiments of Thomson and of Davisson and Germer. They showed that 
electrons (particles) were diffracted like waves with the de Broglie wavelength. 

The de Broglie theory was successful in these instances, but it is incomplete 
and unsatisfying for several reasons. For one, we seldom see particles with a 
unique momentum p\ if the momentum of a particle changes, such as when it is 
acted upon by an external force, its wavelength must change, but the de Broglie 
relationship lacks the capability to enable computation of the dynamical behavior 
of the waves. For this we need a more complete mathematical theory, which was 
supplied by Schrodinger in 1925 and which we review in Section 2 of this 
chapter. A second objection to the de Broglie theory is its reliance on classical 
concepts and terminology. “Particle” and “wave” are mutually exclusive sorts of 
behaviors, but the de Broglie relationship involves classical particles with uniquely 
defined momenta and classical waves with uniquely defined wavelengths. A 
classical particle has a definite position in space. Now, according to de Broglie, 
that localized particle is to be represented by a pure wave that extends throughout 
all space and has no beginning, end, or easily identifiable “position.” 

The solution to this dilemma requires us to discard the classical idea of 
“particle” when we enter the domain of quantum physics. The size of a classical 
particle is the same in every experiment we may do; the “size” of a quantum 
particle varies with the experiment we perform. Quantum physics forces us to 
sacrifice the objective reality of a concept such as “size” and instead to substitute 
an operational definition that depends on the experiment that is being done. Thus 
an electron may have a certain size in one experiment and a very different size in 
another. Only through this coupling of the observing system and the observed 
object can we define observations in quantum physics. A particle, then, is 
localized within some region of space of dimension Ax. It is likely to be found in 
that region and unlikely to be found elsewhere. The dimension Ax of an electron 
is determined by the kind of experiment we do — it may be the dimension of a 
block of material if we are studying electrical conduction in solids, or the 
dimension of a single atom if we are studying atomic physics, or of a nucleus if 
we are studying /? decay. The wave that characterizes the particle has large 
amplitude in the region Ax and small amplitude elsewhere. The single de Broglie 
wave corresponding to the unique momentum component p x had a large ampli- 
tude everywhere; thus a definite momentum (wavelength) corresponds to a 
completely unlocalized particle. To localize the particle, we must add (superpose) 
other wavelengths corresponding to other momenta p x , so that we make the 
resultant wave small outside the region Ax. We improve our knowledge of Ax at 
the expense of our knowledge of p x . The very act of confining the particle to Ax 
destroys the precision of our knowledge of p x and introduces a range of values A p x . 
If we try to make a simultaneous determination of x and p x , our result will show 
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that each is uncertain by the respective amounts Ax and A p x , which are related 
by the Heisenberg uncertainty relationship 

h 

A* A p x >- (2.1) 

with similar expressions for the y and z components. (The symbol h , read as 
“h-bar,” is h/lir where h is Planck’s constant.) The key word here is “simulta- 
neous” — we can indeed measure x with arbitrarily small uncertainty (Ax = 0) if 
we are willing to sacrifice all simultaneous knowledge of the momentum of the 
particle. Having made that determination, we could then make an arbitrarily 
precise measurement of the new momentum (A p x = 0), which would simulta- 
neously destroy our previous precise knowledge of its position. 

We describe the particle by a “wave packet,” a collection of waves, repre- 
senting a range of momenta A p x around p x , with an amplitude that is reasonably 
large only within the region Ax about x. A particle is localized in a region of 
space defined by its wave packet; the wave packet contains all of the available 
information about the particle. Whenever we use the term “particle” we really 
mean “wave packet”; although we often speak of electrons or nucleons as if they 
had an independent existence, our knowledge of them is limited by the uncer- 
tainty relationship to the information contained in the wave packet that describes 
their particular situation. 

These arguments about uncertainty hold for other kinds of measurements as 
well. The energy E of a system is related to the frequency v of its de Broglie wave 
according to E = hv. To determine E precisely, we must observe for a sufficiently 
long time interval A* so that we can determine v precisely. The uncertainty 
relationship in this case is 



h 

AEAt>- (2.2) 

If a system lives for a time At, we cannot determine its energy except to within an 
uncertainty A E. The energy of a system that is absolutely stable against decay 
can be measured with arbitrarily small uncertainty; for all decaying systems there 
is an uncertainty in energy, commonly called the energy “width.” 

A third uncertainty relationship involves the angular momentum. Classically, 
we can determine all three components i x , t z of the angular momentum 
vector i. In quantum mechanics, when we try to improve our knowledge of one 
component, it is at the expense of our knowledge of the other two components. 
Let us choose to measure the z component, and let the location of the projection 
of ( in the xy plane be characterized by the azimuthal angle <p. Then 

K A<J> > \ (2.3) 

If we know t z exactly, then we know nothing at all about <f>. We can think of / as 
rotating or precessing about the z axis, keeping £ z fixed but allowing all possible 
{ x and so that <f> is completely uncertain. 
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2.2 PRINCIPLES OF QUANTUM MECHANICS 



The mathematical aspects of nonrelativistic quantum mechanics are determined 
by solutions to the Schrodinger equation. In one dimension, the time-independent 
Schrodinger equation for a particle of mass m with potential energy V(x) is 

h 2 d 2 xp 

~T~T2 + V(x)*(x) = Et(x) (2.4) 

2 m dx 



where xp(x) is the Schrodinger wave function. The wave function is the mathe- 
matical description of the wave packet. In general, this equation will have 
solutions only for certain values of the energy E ; these values, which usually 
result from applying boundary conditions to xp(x), are known as the energy 
eigenvalues. The complete solution, including the time dependence, is 



^F(x, t) = \p{x) e~ ioit 



(2.5) 



where co = E/h. 

An important condition on the wave function is that xp and its first derivative 
dxp/dx must be continuous across any boundary; in fact, the same situation 
applies to classical waves. Whenever there is a boundary between two media, let 
us say at a: = a, we must have 



and 



lim [xp(a + e) — i p(a — e)] =0 (2.6 a) 

t'—>Q 



lim 

e — »0 



dxp 

dx 



(-) 

\ dx J x = a- e 



(2.6 b) 



It is permitted to violate condition 2.6 b if there is an infinite discontinuity in 
V(x); however, condition 2.6 a must always be true. 

Another condition on xp, which originates from the interpretation of probabil- 
ity density to be discussed below, is that xp must remain finite. Any solution for 
the Schrodinger equation that allows xp to become infinite must be discarded. 

Knowledge of the wave function ^(x, t) for a system enables us to calculate 
many properties of the system. For example, the probability to find the particle 
(the wave packet) between x and x + dx is 

P(x) dx = ^*(x , t) ^(x, t) dx (2.7) 

where x k* is the complex conjugate of T. The quantity 'F* X F is known as the 
probability density. The probability to find the particle between the limits x x and 
x 2 is the integral of all the infinitesimal probabilities: 

P = (2.8) 

x l 

The total probability to find the particle must be 1: 

P ***<& = 1 (2.9) 

- OC 



This condition is known as the normalization condition and in effect it determines 
any multiplicative constants included in 'F. All physically meaningful wave 
functions must be properly normalized. 
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Any function of x, /(x), can be evaluated for this quantum mechanical system. 
The values that we measure for /(x) are determined by the probability density, 
and the average value of /(x) is determined by finding the contribution to the 
average for each value of x: 



(/> = ( 2 . 10 ) 

Average values computed in this way are called quantum mechanical expectation 
values . 

We must be a bit careful how we interpret these expectation values. Quantum 
mechanics deals with statistical outcomes, and many of our calculations are really 
statistical averages. If we prepare a large number of identical systems and 
measure /(x) for each of them, the average of these measurements will be (/>. 
One of the unsatisfying aspects of quantum theory is its inability to predict with 
certainty the outcome of an experiment; all we can do is predict the statistical 
average of a large number of measurements. 

Often we must compute the average values of quantities that are not simple 
functions of x. For example, how can we compute ( p x )l Since p x is not a 
function of x, we cannot use Equation 2.10 for this calculation. The solution to 
this difficulty comes from the mathematics of quantum theory. Corresponding to 
each classical variable, there is a quantum mechanical operator. An operator is a 
symbol that directs us to perform a mathematical operation, such as exp or sin or 
d/dx. We adopt the convention that the operator acts only on the variable or 
function immediately to its right, unless we indicate otherwise by grouping 
functions in parentheses. This convention means that it is very important to 
remember the form of Equation 2.10; the operator is “sandwiched” between 'P* 
and 'P, and operates only on 'P. Two of the most common operators encountered 
in quantum mechanics are the momentum operator, p x = —ihd/dx and the 
energy, E = ihd/dt. Notice that the first term on the left of the Schrodinger 
equation 2.4 is just p x /2m, which we can regard as the kinetic energy operator. 
Notice also that the operator E applied to ^P(x, t) in Equation 2.5 gives the 
number E multiplying ^(x, t). 

We can now evaluate the expectation value of the x component of the 
momentum: 

( Px ) = f^-ih-E^dx 

- -ih /V* — dx (2.11) 

J ox 

One very important feature emerges from these calculations: when we take the 
complex conjugate of 'P as given by Equation 2.5, the time-dependent factor 
becomes e + lui \ and therefore the time dependence cancels from Equations 
2.7-2.11. None of the observable properties of the system depend on the time. 
Such conditions are known for obvious reasons as stationary states ; a system in a 
stationary state stays in that state for all times and all of the dynamical variables 
are constants of the motion. This is of course an idealization — no system lives 
forever, but many systems can be regarded as being in states that are approxi- 
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mately stationary. Thus an atom can make a transition from one “stationary” 
excited state to another “stationary” state. 

Associated with the wave function 'f r is the particle current density j: 





( 2 . 12 ) 



This quantity is analogous to an electric current, in that it gives the number of 
particles per second passing any point x. 

In three dimensions, the form of the Schrodinger equation depends on the 
coordinate system in which we choose to work. In Cartesian coordinates, the 
potential energy is a function of (x, y 9 z) and the Schrodinger equation is 

d 2 xp d 2 \p d 2 ip \ 

+ jyi + + v ( x > Z )’H JC > y> z ) = Eipi*, y, z ) (2.13) 

The complete time-dependent solution is again 

*(*, y, z, t) = Ip(x, y, z) e~ ,ut (2.14) 

The probability density 'F* 'k now gives the probability per unit volume; the 
probability to find the particle in the volume element dv = dx dy dz at x, y, z is 

Pdv = y*<trdv (2.15) 




To find the total probability in some volume V, we must do a triple integral over 
x, y , and z. All of the other properties discussed above for the one-dimensional 
system can easily be extended to the three-dimensional system. 

Since nuclei are approximately spherical, the Cartesian coordinate system is 
not the most appropriate one. Instead, we must work in spherical polar coordi- 
nates (r, #, <f>), which are shown in Figure 2.1. In this case the Schrodinger 
equation is 

h 2 d 2 \p 2 dip 1 d / diM 1 d 2 xp 

2m dr 2 r dr r 2 sin# d6 \ dO j r 2 sin 2 6 d<p 2 
+ V(r, e , C p)t(r 9 e , 4>) = Eip(r , fi, <j>) (2.16) 




Figure 2.1 Spherical polar coordinate system, showing the relationship to Carte 
sian coordinates. 
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All of the previous considerations apply in this case as well, with the volume 
element 

do = r 2 sin6 dr d8d<l> (2.17) 

The following two sections illustrate the application of these principles, first 
with the mathematically simpler one-dimensional problems and then with the 
more physical three-dimensional problems. 

2.3 PROBLEMS IN ONE DIMENSION 



The Free Particle 



For this case, no forces act and we take V(x) = 0 everywhere. We can then 
rewrite Equation 2.4 as 



d 2 \p 

dx 2 



2 mE 




(2.18) 



The solution to this differential equation can be written 

^(x) = A ' sin kx + B f cos kx (2.19) 

or, equivalently 

$(x) =Ae ikx + Be~ ikx (2.20) 

where k 2 = 2 mE/h 2 and where A and B (or A’ and B') are constants. 

The time-dependent wave function is 

*(*, t)=A e iikx ~ at) + Be~ i(kx+Ut) (2.21) 

The first term represents a wave traveling in the positive * direction, while the 
second term represents a wave traveling in the negative x direction. The intensi- 
ties of these waves are given by the squares of the respective amplitudes, \A\ 2 
and \B\ 2 . Since there are no boundary conditions, there are no restrictions on the 
energy E\ all values of E give solutions to the equation. The normalization 
condition 2.9 cannot be applied in this case, because integrals of sin 2 or cos 2 do 
not converge in x = - oo to 4- oo. Instead, we use a different normalization 
system for such constant potentials. Suppose we have a source such as an 
accelerator located at x = — oo, emitting particles at a rate I particles per 
second, with momentum p = hk in the positive * direction. Since the particles 
are traveling in the positive .x direction, we can set B to zero — the intensity of 
the wave representing particles traveling in the negative * direction must vanish 
if there are no particles traveling in that direction. The particle current is, 
according to Equation 2.12, 

hk , / x 

j=—\A | 2 (2.22) 

m 

which must be equal to the current of I particles per second emitted by the 
source. Thus we can take A = yjml/hk . 
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Step Potential, E > V 0 



The potential is 



V(x) = 0 x<0 

= V 0 x > 0 



(2.23) 



where V 0 > 0. Let us call jc < 0 region 1 and * > 0 region 2. Then in region 1, 
the Schrodinger equation is identical with Equation 2.18 and the solutions are 
given by Equation 2.20 with k = k x = ]/2mE/h 2 . In region 2, the Schrodinger 
equation is 



d 2 ^ 2 2 m(E - V 0 ) 

dx 2 ~ h 2 ^ 


(2.24) 


Since £ > V 0 , we can write the solution as 




i// 2 = CV* 2 * + De~ ik '- X 


(2.25) 


where k 2 = )j2m(E - V Q )/h 2 . 

Applying the boundary conditions at jc = 0 gives 




A+B=C+D 


( 2 . 26 a) 


from Equation 2.6a, and 




k x {A - B) = k 2 (C - D) 


( 2 . 26 b) 



from Equation 2.6b. 

Let’s assume that particles are incident on the step from a source at jc = - oo. 
Then the A term in represents the incident wave (the wave in jc < 0 traveling 
toward the step at jc = 0), the B term in \p x represents the reflected wave (the 
wave in jc < 0 traveling back toward jc = - oo), and the C term in \p 2 represents 
the transmitted wave (the wave in jc > 0 traveling away from jc = 0). The D term 
cannot represent any part of this problem because there is no way for a wave to 
be moving toward the origin in region 2, and so we eliminate this term by setting 
D to zero. Solving Equations 2.26a and 2.26b , we have 

1 - k 2 /k x 

B = A — 1 

1 + ^2/^1 




C = 



2 

A 

1 + k 2 /k x 



(2.28) 



The reflection coefficient R is defined as the current in the reflected wave divided 
by the incident current: 



R 



J reflected 



-/incident 



and using Equation 2.22 we find 

*=^.(1^1) 

\A\ 2 \ 1 + k 2 /k 1 J 



(2.29) 



(2.30) 



The transmission coefficient T is similarly defined as the fraction of the incident 
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Figure 2.2 The wave function of a particle of energy £ encountering a step of 
height V 0 for the case E > V 0 . The de Broglie wavelength changes from to \ 2 
when the particle crosses the step, but ^ and d-ty/dx are continuous at x - 0. 



current that is transmitted past the boundary: 



and thus 






^transmitted 
J incident 



ki\C^ _ 4fe 2 /fc 1 
k y \A\ 2 (l + k 2 /k y ) 2 



(2.31) 



(2.32) 



Notice that R + T = 1, as expected. The resulting solution is illustrated in 
Figure 2.2. 

This is a simple example of a scattering problem. In Chapter 4 we show how 
these concepts can be extended to three dimensions and applied to nucleon- 
nucleon scattering problems. 



Step Potential, E < V Q 

In this case, the potential is still given by Equation 2.23, and the solution for 
region 1 (x < 0) is identical with the previous calculation. In region 2, the 
Schrodinger equation gives 

d 2 \bj 2m 

(2.33) 

which has the solution 

ip 2 = Ce klX + De~ k * x (2.34) 

where k 2 = ^2 m(V 0 - E)/h 2 . Note that for constant potentials, the solutions 
are either oscillatory like Equation 2.19 or 2.20 when E > V 0 , or exponential like 
Equation 2.34 when E < V 0 . Although the mathematical forms may be different 
for nonconstant potentials V(x), the general behavior is preserved: oscillatory 
(though not necessarily sinusoidal) when E > V(x) and exponential when E < 
V(x). 

This solution, Equation 2.34, must be valid for the entire range x > 0. Since 
the first term would become infinite for x -> oo, we must set C = 0 to keep the 
wave function finite. The D term in ip 2 illustrates an important difference 
between classical and quantum physics, the penetration of the wave function into 
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V V V v 0 

J i 

*=o 

Figure 2.3 The wave function of a particle of energy E encountering a step of 
height V 0 , for the case E < V 0 . The wave function decreases exponentially in the 
classically forbidden region, where the classical kinetic energy would be negative. 
At x = 0, and d\p/dx are continuous. 

the classically forbidden region . All (classical) particles are reflected at the 
boundary; the quantum mechanical wave packet , on the other hand, can penetrate 
a short distance into the forbidden region. The (classical) particle is never directly 
observed in that region; since E < V 0 , the kinetic energy would be negative in 
region 2. The solution is illustrated in Figure 2.3 



(2.35) 



Barrier Potential, E > V Q 

The potential is 

V{x) = 0 jc < 0 

= V 0 0 < x < a 
= 0 x > a 

In the three regions 1, 2, and 3, the solutions are 

i/q = Ae ik ' x + Be~ ik ' x 
= Ce lk2X + De~ ik * x 
= Fe ik > x + Ge~ ik i x 

where = k 3 = ^jlmE/h 1 and k 2 = /2 m(E - V Q )/h 2 . 

Using the continuity conditions at jc = 0 and at jc = a, and assuming again 
that particles are incident from jc = — oo (so that G can be set to zero), after 
considerable algebraic manipulation we can find the transmission coefficient 
T= \F\ 2 /\A\ 2 -. 

1 

T = : — (2.37) 



(2.36) 



1 

1 + - 



V 2 



4 E(E 

The solution is illustrated in Figure 2.4. 



Vo) 



sin 2 k 7 a 



Barrier Potential, E < V 0 

For this case, the ip l and ip 3 solutions are as above, but \p 2 becomes 

= Ce k ’ x + De~ k ' x (2.38) 



where now k 2 = \2m ( V 0 - E)/h 2 . Because region 2 extends only from x = 0 
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Figure 2.4 The wave function of a particle of energy E > V Q encountering a 
barrier potential. The particle is incident from the left. The wave undergoes reflec- 
tions at both boundaries, and the transmitted wave emerges with smaller amplitude. 



to x = a, the question of an exponential solution going to infinity does not arise, 
so we cannot set C or D to zero. 

Again, applying the boundary conditions at x = 0 and x = a permits the 
solution for the transmission coefficient: 



T = 



1 



1 + 



1 v o 

4 E(V 0 -E) 



sinh 2 k 2 a 



(2.39) 



Classically, we would expect T = 0— the particle is not permitted to enter the 
forbidden region where it would have negative kinetic energy. The quantum wave 
can penetrate the barrier and give a nonzero probability to find the particle 
beyond the barrier. The solution is illustrated in Figure 2.5. 

This phenomenon of barrier penetration or quantum mechanical tunneling has 
important applications in nuclear physics, especially in the theory of a decay, 
which we discuss in Chapter 8. 



s\ A /> 

E ■ / \ ; 






E \ j % f 

v \J 







x = 0 x = a 



Figure 2.5 The wave function of a particle of energy E < V 0 encountering a 
barrier potential (the particle would be incident from the left in the figure). The 
wavelength is the same on both sides of the barrier, but the amplitude beyond the 
barrier is much less than the original amplitude. The particle can never be ob- 
served, inside the barrier (where it would have negative kinetic energy) but it can 
be observed beyond the barrier. 
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The Infinite Well 



The potential is (see Figure 2.6) 

V(x) = oo x < 0, x > a 

= 0 0 < x < a (2.40) 



That is, the particle is trapped between x = 0 and x — a. The walls at * = 0 and 
.x = a are absolutely impenetrable; thus the particle is never outside the well and 
xp = 0 for jc < 0 and for jc > a. Inside the well, the Schrodinger equation has the 
form of Equation 2.18, and we will choose a solution in the form of Equation 
2.19: 



xp = A sin kx + B cos kx 



(2.41) 



The continuity condition on \p at x = 0 gives \p(0) = 0, which is true only for 
B = 0. At jc = a, the continuity condition on xp gives 

A sin ka = 0 (2.42) 



The solution A = 0 is not acceptable, for that would give ip = 0 everywhere. 
Thus sin ka = 0, or 



and 



ka = mr n — 1,2,3,... 

h 2 k 2 h 2 n t 2 „ 

p = = n 2 

n 2m Ima 2 



(2.43) 

(2.44) 



To OO 



To OO 




Figure 2.6 A particle moves freely in the one-dimensional region 0 < x < a but is 
excluded completely from x < 0 and x > a. A bead sliding without friction on a wire 
and bouncing elastically from the walls is a simple physical example. 
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Here the energy is quantized — only certain values of the energy are permitted. 
The energy spectrum is illustrated in Figure 2.7. These states are bound states , in 
which the potential confines the particle to a certain region of space. 

The corresponding wave functions are 



tn = 




(2.45) 




Excited 

states 




E = 16E 0 





Figure 2.7 The permitted energy levels of the one-dimensional infinite square 
well. The wave function for each level is shown by the solid curve, and the shaded 
region gives the probability density for each level. The energy E 0 is h 2 * 2 / 2ma 2 . 
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where the constant A has been evaluated using Equation 2.9. The probability 
densities | ^ | 2 of some of the lower states are illustrated in Figure 2.7. 

The Finite Potential Well 

For this case we assume the well has depth V 0 between + a/2 and —a/2: 

V(x) = V 0 \x\ > a/2 

= 0 1*1 < a/2 (2.46) 

We look for bound-state solutions, with E < V 0 . The solutions are 
^ = Ae k ' x + Be~ k ' x x < -a/2 

= Ce iklX + De~ ik * x — a/2 < x < a/2 (2.47) 

if> 3 = Fe k ' x + Ge~ k ' x x > a/2 

where k x = y2 m(V {) — E)/h 2 and k 2 = / 2mE/h 2 . To keep the wave function 
finite in region 1 when jc -> — oo, we must have B — 0, and to keep it finite in 
region 3 for x -> 4- oo, we require F = 0. 

Applying the continuity conditions at x = —a/2 and at jc = +a/ 2, we find 
the following two relationships: 

k 2 ci 

k 2 tan — = k x (2.48a) 

or 

k 7 a 

-k 2 cot — = k x (2.486) 

These transcendental equations cannot be solved directly. They can be solved 
numerically on a computer, or graphically. The graphical solutions are easiest if 
we rewrite Equations 2.48 in the following form: 

a tan a = (P 2 — a 2 ) 1/2 (2.49a) 

- acota = (P 1 - a 2 ) 1/2 (2.496) 

where a = k 2 a/ 2 and P = (mV 0 a 2 /2h 2 )' The right side of these equations 

defines a circle of radius P, while the left side gives a tangentlike function with 
several discrete branches. The solutions are determined by the points where the 
circle intersects the tangent function, as shown in Figure 2.8. Therefore , the 
number of solutions is determined by the radius P, and thus by the depth V 0 of the 
well. (Contrast this with the infinite well, which had an infinite number of bound 
states.) For example, when P < 77/2, there is only one bound state. For m/2 < 
P < m there are two bound states. Conversely, if we studied a system of this sort 
and found only one bound state, we could deduce some limits on the depth of the 
well. As we will discuss in Chapter 4, a similar technique allows us to estimate the 
depth of the nuclear potential, because the deuteron, the simplest two-nucleon 
bound system, has only one bound state. 
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a 

(a) 




( 6 ) 

Figure 2.8 (a) The graphical solution of Equations 2.49a and 2.496. For the case 
of P = 6 (chosen arbitrarily) there are four solutions at a = 1 .345, 2.679, 3.985, and 
5.226. ( b ) The wave functions and probability densities (shaded) for the four states. 
(Compare with the infinite well shown in Figure 2.7.) 



The Simple Harmonic Oscillator 

Any reasonably well-behaved potential can be expanded in a Taylor series about 
the point x 0 : 



V(x) = V(x 0 ) + 



IdV 

\ ^7 



( X - x 0 ) + 

JC = JC„ 



1 

2 



d 2 V 

dx 2 



(x- x 0 ) 2 + ••• 

(2.50) 



If x 0 is a potential minimum, the second term in the series vanishes, and since 
the first term contributes only a constant to the energy, the interesting term is the 
third term. Thus to a first approximation, near its minimum the system behaves 
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Table 2.1 Sample Wave Functions of the One-Dimensional 
Simple Harmonic Oscillator 



n 


E„ 


'f'Jx) 


0 




v -l/4 e -a 2 x 2 /2 


1 


§/j<o 0 


2~ I/ 2 ,," >/ 4 (2«x) e -“ 2 * 2 /2 


2 


§/i<Oo 


2 -3/2 ff -l/4 ( 4 a 2 x 2 _ 2) e -« 2 * 2 /2 


3 


jhu 0 


(l/4v / 3ff 1 / 4 )(8« 3 x 3 - Vlax) e ~ a2x2/1 


4 




(I/ 81/6 77 1/4 )(16a 4 jc 4 - 48a 2 x 2 + 12) e ”“ 2j[2 / 2 




E„ = 


hu 0 (n + j) 




'PAx) = 


(2"n\^y 1/2 H n (ax) e~ a2x2/2 
where H n (ax) is a Hermite polynomial 



like a simple harmonic oscillator, which has the similar potential \k(x - x 0 ) 2 . 
The study of the simple harmonic oscillator therefore is important for under- 
standing a variety of systems. 

For our system, we choose the potential energy 

V(x) = \kx 2 (2.51) 

for all x. The Schrodinger equation for this case is solved through the substitu- 
tion t//(x) = h(x) e~ a2x2/ 2 , where a 2 = {km/h. The function h(x) turns out to 
be a simple polynomial in x. The degree of the polynomial (the highest power of 
x that appears) is determined by the quantum number n that labels the energy 
states, which are also found from the solution to the Schrodinger equation: 

E n = ho) 0 (n + \) n = 0, 1,2,3,... (2.52) 

where co 0 = {k/m , the classical angular frequency of the oscillator. Some of the 
resulting wave functions are listed in Table 2.1, and the corresponding energy 
levels and probability densities are illustrated in Figure 2.9. Notice that the 
probabilities resemble those of Figure 2.8; where E > F, the solution oscillates 
somewhat sinusoidally, while for E < V (beyond the classical turning points 
where the oscillator comes to rest and reverses its motion) the solution decays 
exponentially. This solution also shows penetration of the probability density 
into the classically forbidden region. 

A noteworthy feature of this solution is that the energy levels are equally 
spaced. Also, because the potential is infinitely deep, there are infinitely many 
bound states. 

Summary 

By studying these one-dimensional problems, we learn several important details 
about the wave properties of particles. 

1. Quantum waves can undergo reflection and transmission when they encoun- 
ter a potential barrier; this behavior is very similar to that of classical waves. 

2. A wave packet can penetrate into the classically forbidden region and appear 
beyond a potential barrier that it does not have enough energy to overcome. 
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Figure 2-9 The lowest few energy levels and corresponding probability densities 
of the harmonic oscillator. 



3. Wave functions oscillate whenever E > V(x) and decay exponentially 
whenever E < V{x). 

4. When a potential confines a particle to a region of space, bound-state wave 
functions can result. The particle is permitted only a set of discrete energy 
values; the number of allowed energy values is determined by the depth of 
the potential well. 



2.4 PROBLEMS IN THREE DIMENSIONS 



The Infinite Cartesian Well 

We begin with a problem in Cartesian coordinates that illustrates an important 
feature present in three-dimensional problems but not in one-dimensional prob- 
lems. The potential in this case is 

V(x,y, z) = 0 0 <x<a, 0 <y<a, 0 <z<a (2 53) 

= oo x < 0, x > a, y < 0, y > a, z < 0, z > a 

The particle is thus confined to a cubical box of dimension a . Beyond the 
impenetrable walls of the well, \p = 0 as before. Inside the well, the Schrodinger 
equation is 



h 



d 2 \p d 2 \p d 2 \p 

2m \ dx 2 dy 2 dz 2 



= E\P(x,y,z) (2.54) 

The usual procedure for solving partial differential equations is to try to find a 
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(4,2,1) 






(3,3,1) 


19 




(4,1,1) 






(3,2,2) 


17 




(3,2,1) 


14 




(2,2,2) 


12 




(3,1,1) 




3 


(2,2,1) 






(2,1,1) 


6 


3 ni " ,,n " 


(1,1,1) 
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Degeneracy {n x , n y , n z ) Energy 





Figure 2.10 Energy levels of a particle confined to a three-dimensional cubical 
box. The energy is given in units of E 0 = /Fir 2 / 2ma 2 . 



separable solution, with \p(x, y, z) = X(x)Y(y) Z(z), where X, Y, and Z are 
each functions of a single variable. We will skip the mathematical details and give 
only the result of the calculation: 



y, z) 




n y 7ry n y mz 

sin — sin — — 

a a 



(2.55) 



2 _ 2 



n x n v n : 



h 2 *TT 
2 ma 



:(«* + n l + n l ) 



(2.56) 



where n x> n v , and n z are independent integers greater than zero. The lowest 
state, the ground state, has quantum numbers ( n x , n v , n z ) = (1,1,1). Its prob- 
ability distribution would show a maximum at the center of the box (x = y = 
z = a/2), falling gradually to zero at the walls like sin 2 . 

The first excited state has three possible sets of quantum numbers: (2,1,1), 
(1, 2, 1), and (1, 1, 2). Each of these distinct and independent states has a different 
wave function, and therefore a different probability density and different expecta- 
tion values of the physical observables, but they all have the same energy. This 
situation is known as degeneracy ; the first excited state is threefold degenerate. 
Degeneracy is extremely important for atomic structure since it tells us how 
many electrons can be put in each atomic subshell. We will soon discuss its 
similar role in the nuclear shell model. 

Figure 2.10 shows the lower portion of the energy spectrum of excited states. 
Notice that the spacing and ordering do not have the regularity of the one-dimen- 
sional problem. 
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The Infinite Spherical Well 



If we work in spherical coordinates with a potential that depends only on r (not 
on 0 or <J>), another new feature arises that will be important in our subsequent 
investigations of nuclear structure. When we search for separable solutions, of the 
form \f/(r, 0, <J>) = R(r) 0(0) 0(<J>), the central potential V(r) appears only in 
the radial part of the separation equation, and the angular parts can be solved 
directly. The differential equation for $(<£) is 

d 2 0 

-^T + m) * = 0 (2.57) 

where is the separation constant. 

The solution is 



1 d 
sin# dO 



= 




. . . The equation for 0(0) is 


/ d@\ 


m 2 




1 - -T7 
sin 0 



0 = 0 



(2.58) 



(2.59) 



where £ = 0, 1, 2, 3, . . . and m tf = 0, ±1, ±2, . . . , ±£. The solution 0^(0) can 
be expressed as a polynomial of degree £ in sin 0 or cos 0. Together, and 
normalized, 0^(0) and 0^(0) give the spherical harmonics Y emt (0, <#>), some 
examples of which are listed in Table 2.2. These functions give the angular part 
of the solution to the Schrodinger equation for any central potential V(r). For 
example, it is these angular functions that give the spatial properties of atomic 
orbitals that are responsible for molecular bonds. 

For each potential F(r), all we need to do is to find a solution of the radial 
equation 



2m 



d 2 R 2 dR 

~Z~ + 

dr r dr 



V{r) + 



1 )n 2 

2 mr 2 



R = ER 



(2.60) 



Table 2.2 Spherical Harmonics for Some Low £ Values 



£ 


m e 




0 


0 


( l/4w) 1/2 


1 


0 


(3/4tt) 1/2 cos 0 


1 


+ 1 


+ (3/87r) 1/2 sin 0 e ±i ' 1 ’ 


2 


0 


(S/16ir) 1/2 (3cos 2 0 - 1) 


2 


±1 


: F(15/87r) 1/2 sin0cos0 e*'* 


2 


+ 2 


(15/327r) 1 / 2 sin 2 0 e ±2i * 



*»,(*) = 









fa 
2£ + 1 (£ — 



1/2 



P?'(0) 



2 (£+ m^)\ 

where P" l '{0) is the associated Legendre polynomial 
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Table 2-3 Spherical Bessel Functions — 
Sample Expressions and Limits 



Jo( kr ) = 
Mkr) = 

Jiikr) = 

M kr ) s 

jf(kr) s 



sin kr 



kr 

sin kr cos kr 



(kr) kr 

3 sin kr 3 cos kr sin kr 



(kr) 3 (kr) 2 

(krY 

1 • 3 • 5 • • • (2/ + 1) 
sin( kr — *?V/2) 



kr 

kr^O 



kr 

Si 



kr oo 



■"<*'> ‘("l) (7^) AW 



The *f(/+ 1) term is generally written as an addition to the potential; it is called 
the “centrifugal potential” and it acts like a potential that keeps the particle away 
from the origin when ( > 0. 

As an example, we consider the case of the infinite spherical well, 

V(r) = 0 r<a 

= 00 r > a (2.61) 

We require again that R(r ) = 0 for r > a, since the walls of the infinite well are 
impenetrable. Inside the well, the solution to Equation 2.60 for V = 0 can be 
expressed in terms of the oscillatory functions known as spherical Bessel functions 
js(kr), some of which are listed in Table 2.3. To find the energy eigenvalues, we 
proceed exactly as in the one-dimensional problem and apply the continuity 
condition on ^ at r — a. This gives 

jAka) = 0 (2.62) 

This is in effect a transcendental equation, which must be solved numerically. 
Tables of the spherical Bessel functions are available that can be consulted to 
find the zeros for any given value of *f.* For example, we consider the case S= 0. 
From the tables we find y 0 (x) = 0 at x = 3.14, 6.28, 9.42, 12.57, and so on. For 
tf = 1, the first few zeros of j\(x) are at x = 4.49, 7.73, 10.90, 14.07. Since 
E = h 2 k 2 /lm, we can then solve for the allowed values of the energies. Repeat- 
ing this process for 2, f— 3, and so on, we would be able to construct a 
spectrum of the energy states, as is shown in Figure 2.11. As in the case of the 
Cartesian well, the regularity of the one-dimensional problem is not present. 
Also, note that the levels are again degenerate — since the energy depends only on 
i , the wave functions with different m e values all have the same energy. Thus, in 
the case of the level with ^=2, the possible wave functions are j 2 (kr)Y 22 (0, <f>), 



M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions (New York: Dover, 1965). 
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Figure 2.1 1 Energy levels of a particle confined to a three-dimensional spherical 
container. The energy is given in units of E 0 = h 2 ^ 2 /2ma 2 . Compare with the 
spacings and degeneracies of Figure 2.10. The quantum number n does not arise 
directly in the solution in this case; it serves to number the states of a given 



j 2 (kr)Y 2l (6, <f>), j 2 (kr )Y 20 (6, $), j 2 (kr)Y 2 _ 1 (8, <f>), and j 2 (kr)Y 2 _ 2 (8, <f>), for a 
fivefold degeneracy. In fact, since m e is restricted to the values 0, ±1, 
± 2, . . . , ± rf, there are exactly 2/+ 1 possible Y*m e for a given /, and thus each 
level has a degeneracy of 2/+ 1. (This situation is very similar to the case of 
electronic orbits in atoms, in which there is also a central potential. The capacity 
of each atomic subshell contains the factor of 2/+ 1, which likewise arises from 
the nif degeneracy.) 

The probability to locate the particle in a volume du is given by \\p\ 2 dv, where 
the volume element was given in Equation 2.17. Such three-dimensional distribu- 
tions are difficult to represent graphically, so we often consider the radial and 
angular parts separately. To find the radial probability density, which gives the 
probability to find the particle between r and r + dr averaged over all angles, we 
integrate the probability density over 6 and (j>: 

P(r)dr= (W 2 dv 

= r 2 \R(r)\ 2 dr J sin 6 dO j d<f>\Y^ m f (2.63) 

The spherical harmonics Y tmt are themselves normalized, so that the integral 
gives 1, and thus 

P(r) = r 2 |tf(r)| 2 (2.64) 

Some sample radial probability distributions for the infinite well are shown in 
Figure 2.12. 

The angular dependence of the probability density for any central potential is 
given by | ( 6 , </>) | 2 , some samples of which are illustrated in Figure 2.13. 
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Figure 2.12 The left side shows, for some of the lower energy levels, the 
unnormalized y>(/cr), adjusted so that j^ka) = 0. The right side shows the corre- 
sponding, normalized radial probability density, r 2 R 2 . Note that all j vanish at the 
origin except j 0 , and that all probability densities vanish at r = 0. Also, note how the 
“centrifugal repulsion” pushes corresponding maxima in P(r) away from the origin 
as t increases. 




Figure 2.13 Spatial probability distributions resulting from the Y, m . The three- 
dimensional representations can be obtained by rotating each figure" about the z 
a> 
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Table 2.4 Sample Radial Wave Functions for Three-Dimensional 
Simple Harmonic Oscillator 



n 


£ 


E„ 


R(r) 


0 


0 


o 


(2« 3 /2 /w ‘/ 4 ) e -«V/2 


1 


1 


t&«o 


(2 o? /2 yfl //3 tt 1 / 4 )(ar) e - " 2 ' 2 / 2 


2 


0 


1^0 


(2a 3/2 v^/V3V /4 )(§ - a 2 r 2 ) e~ aV/2 


2 


2 


O 

3 

•« 


(4a 3 / 2 //L5 ir I/4 )(a 2 r 2 ) e~ a2r2/2 


3 


1 


f*«l 0 


(4a 3/2 //T5 w 1/4 )(far - a 3 r 3 ) e ~“ v /2 


3 


3 


!»«b 


(4a 3/2 v^//i05 w 1/4 )(a 3 r 3 ) e -“ 2 '' 2 /2 


4 


0 




(4a 3/2 v^2//l5w 1 ' /4 )(^ - fa 2 r 2 + fa 4 /- 4 ) <T aV/2 


4 


2 




(4a 3/2 v/2/\/l05w 1/4 )(fa 2 r 2 - a 4 /- 4 ) e -“ 2 '' 2 / 2 


4 


4 


0 


(8a 3/,2 /3v , 105 77 1/4 ) a 4 /- 4 e - ° 2 '' 2 / 2 



Note the similarity in form (polynomial X exponential) between these solutions and those of the 
one-dimensional problem shown in Table 2.1. In this case the polynomials are called Laguerre 
polynomials. The solutions are discussed in J. L. Powell and B. Crasemann, Quantum Mechanics 
(Reading, MA: Addison-Wesley, 1961), Chapter 7. 

The Simple Harmonic Oscillator 

We consider a central oscillator potential, V(r) = \kr 2 . The angular part of the 
solution to the Schrodinger equation is Y, mt X 6 , <j>) for all central potentials, so all 
we need to consider here is the solution to the radial equation. As in the 
one-dimensional case, the solution can be expressed as the product of an 
exponential and a finite polynomial. Some representative solutions are listed in 
Table 2.4, and the corresponding radial probability densities are illustrated in 
Figure 2.14. The general properties of the one-dimensional solutions are also 
present in this case: oscillation in the classically allowed region, and exponential 
decay in the classically forbidden region. 

The energy levels are given by 

E„ = hu 0 (n+§) (2.65) 

where n = 0, 1, 2, 3, . . . . The energy does not depend on £ , but not all £ values 
are permitted. From the mathematical solution of the radial equation, the 
restrictions on £ are as follows: £ can be at most equal to n and takes only even 
or only odd values as n is even or odd. For n = 5, the permitted values of £ are 
1, 3, and 5; for n = 4, the values of £ are 0, 2, and 4. Since the energies do not 
depend on either, there is an additional degeneracy of 2£ + 1 for each £ 
value. Thus the n = 5 level has a degeneracy of [(2 X 1 + 1) + (2 X 3 + 1) + 
(2x5 + 1)] = 21, while the n = 4 level has a degeneracy of [(2x0 + 1) + 
(2 X 2 + 1) + (2 X 4 + 1)] = 15. Figure 2.15 shows some of the energy levels 
and their degeneracies, which are equal to \{n + l)(n + 2). 

The Coulomb Potential 

The attractive Coulomb potential energy also has a simple central form, V(r) = 
— Ze 2 /4'rre 0 r, for the interaction between electrical charges + Ze and —e, such 
as in the case of a single electron in an atom of atomic number Z. The angular 
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Figure 2-14 Radial probabililty densities for some states of the three-dimensional 
harmonic oscillator. The vertical arrows mark the classical turning points. As in 
Figure 2.12, note that P(r) vanishes for r = 0 (but note from Table 2.4 that R(r) is 
nonvanishing at r = 0 only for £= 0). Also note the “centrifugal repulsion” for the 
larger £ values. 



part of the wave function is again given by 8 , <#>), and some of the radial 

wave functions R(r ) are listed in Table 2.5. The energy levels are E n = 
( — mZ 2 e 4 /32'7T 2 c 0 2 h 2 n 2 ) and are shown in Figure 2.16; the radial probability 
density is plotted for several states in Figure 2.17. The relationship between n 
and £ is different from what it was for the oscillator potential: for each level n , 
the permitted values of £ are 0,1,2,... ,(« - 1). The total degeneracy of each 
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Figure 2.15 Some of the lower energy levels of a particle in a central three- 
dimensional oscillator potential. 



Table 2.5 Coulomb (Hydrogenic) Radial Wave Functions 



n 


t 


R(r) 


1 


0 


2(Z/a 0 ) 3/2 e ~ Zr/a ° 


2- 


0 


(Z/2a 0 ) 3/2 (2 - Zr/a 0 ) e- z ' /2 "° 


2 


1 


3“ 1 / 2 (Z/2a 0 ) 3/2 (Zr/a 0 ) e~ Zr / 2fl ° 


3 


0 


f(Z/3a 0 ) 3/2 [3 - 2 Zr/a 0 + 2(Zr/la 0 ) 2 ] e~ Zr ^ 


3 


1 


(4/2 /9)( Z/3a 0 ) 3/2 ( Zr/ a 0 )(l - Zr/6a 0 ) e - z '/ 3 ‘° 


3 


2 


(2/2 /27/F)( Z/3a 0 ) 3/2 (Zr/a 0 ) 2 e - z '/ 3 “» 



The radial wave functions have the mathematical form of associated Laguerre polynomials multiplied 
by exponentials. The Bohr radius a 0 is 4iTe 0 h 2 /me 2 . For a discussion of these solutions, see L. 
Pauling and E. B. Wilson, Introduction to Quantum Mechanics (New York: McGraw-Hill, 1935), 
Chapter 5. 



energy level, including the various ( values and the 2^+1 degeneracy of each, is 



Summary 

The mathematical techniques of finding and using solutions to the Schrodinger 
equation in three dimensions are similar to the techniques we illustrated previ- 
ously for one-dimensional problems. There are two important new features in the 
three-dimensional calculations that do not arise in the one-dimensional calcula- 
tions: (1) The energy levels are degenerate — several different wave functions can 
have the same energy. The degeneracies will have the same effect in the nuclear 
shell model that the i and m e degeneracies of the energy levels of the Coulomb 
potential have in the atomic shell model — they tell us how many particles can 
occupy each energy level. (2) When the potential depends only on r (not on 0 or 
<j>), the wave functions can be assigned a definite angular momentum quantum 
number i . These new features will have important consequences when we discuss 
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r (nm) 




Figure 2.16 The lower energy levels in a Coulomb potential, shown for Z = 1 
(hydrogen atom). The states are labeled with ( n , /); the degeneracies are indicated 
on the left and the energy values on the right. 



nuclear models in Chapter 5. The behavior of angular momentum in quantum 
theory is discussed in the next section. 



2.5 QUANTUM THEORY OF ANGULAR MOMENTUM 

In solutions of the Schrodinger equation for three-dimensional problems, the 
quantum number / plays a prominent role. In atomic physics, for example, it 
serves to label different electron wave functions and to tell us something about 
the spatial behavior of the wave functions. This angular momentum quantum 
number has the same function in all three-dimensional problems involving central 
potentials, where V = V(r). 

In classical physics, the angular momentum f of a particle moving with linear 
momentum p at a location r from a reference point is defined as 

S=rXp (2.66) 

In quantum mechanics, we can evaluate the expectation value of the angular 
momentum by analogy with Equation 2.10. We first consider the magnitude of 
the angular momentum, and for this purpose it is simplest to calculate ( 2 . We 
must first find a quantum mechanical operator for / 2 , as we discussed in Section 
2.2. This can be done simply by replacing the components of p with their 
operator equivalents: p x = —ih d/dx , p y = —ih d/dy , p z = —ih d/dz. Evaluat- 
ing the cross product then gives terms of the form t x = yp z - zp y , and finally 
computing 2 ) = (t 2 + t y 4- t 2s ) gives the remarkably simple result, which is 
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Figure 2-17 Radial probability distributions for a particle in a Coulomb potential 
(hydrogenic atom). The probability vanishes at r = 0, but as before the *?= 0 wave 
functions do not. This property becomes especially important for phenomena that 
depend on the overlap of atomic wave functions with the nucleus— only /= 0 
states contribute substantially to such phenomena (electron capture, hyperfine 
structure, etc.). Why doesn’t the “centrifugal repulsion” appear to occur in this 
case? 



independent of the form of R(r ), 

<Y 2 > = 1) (2.67) 

That is, whenever we have a central potential, which gives a wave function 
R(r)Yf m ^(6, </>), the magnitude of the angular momentum is fixed at the value 
given by Equation 2.67; the angular momentum is a constant of the motion (as it is 
in classical physics for central potentials). The atomic substates with a given i 
value are labeled using spectroscopic notation ; we use the same spectroscopic 
notation in nuclear physics: s for 0, p for 1, and so on. These are 
summarized in Table 2.6. 

When we now try to find the direction of /, we run into a barrier imposed by 
the uncertainty principle: quantum mechanics permits us to know exactly only 

Table 2.6 Spectroscopic Notation 

€ value 0 1 2 3 4 5 6 

Symbol s p d f g h i 
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Figure 2.18 The vector precesses rapidly about the z axis, so that 4 stays 
constant, but 4 and 4 are variable. 

one component of / at a time. Once we determine the value of one component, 
the other two components are completely indeterminate. (This is a fundamental 
limitation, and no amount of trickery can get us around it. It is the very act of 
measuring one component that makes the other two indeterminate. When we 
measure 4 , we force tf y and 4 into indeterminacy; when we then measure 
for the same system, our previous knowledge of t x is destroyed as 4 is now 
forced into indeterminacy.) By convention, we usually choose the z component of 
/ to be determined, and computing (4) as described above, 

«*> = hm e ( 2 . 68 ) 

where m e = 0, +1, + 2, . . . , + (. Notice that |(/ z ) | < \t\ = h\[7J/ + 1) — the z 
component of the vector is always less than its length. If |(4)| = Kl were 
permitted, then we would have exact knowledge of all three components of t (€ x 
and t y would be zero if / were permitted to align with the z axis). The 
conventional vector representation of this indeterminacy is shown in Figure 2.18 
— t rotates or precesses about the z axis keeping ( z fixed but varying £ x and ( , 
The complete description of an electronic state in an atom requires the 
introduction of a new quantum number, the intrinsic angular momentum or spin . 
For the electron, the spin quantum number is s = The spin can be treated as 
an angular momentum (although it cannot be represented in terms of classical 
variables, because it has no classical analog). Thus 

(s 2 ) = h 2 s(s + 1) (2.69) 

{s z ) = hm s (m,= ±£) (2.70) 

It is often useful to imagine the spin as a vector s with possible z components 

+ 

Nucleons, like electrons, have spin quantum numbers of A nucleon moving 
in a central potential with orbital angular momentum f and spin s has a total 
angular momentum 

j = S+s (2.71) 

The total angular momentum j behaves in a manner similar to ( and s: 

( j 2 > = h 2 j(j + !) 

U> = <4 + O = hm j 



(2.72) 

(2.73) 
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Figure 2.19 The coupling of orbital angular momentum / to spin angular 
momentum s giving total angular momentum j. (a) Coupling giving j = S + The 
vectors i and s have definite lengths, as does j. The combined / and s vectors 
rotate or precess about the direction of y; in this coupling the z components 4 and 
s z thus do not have definite values. The vector y precesses about the z direction so 
that j z has a definite value. ( b ) The similar case of j = In interpreting both 

figures, keep in mind that all such representations of vectors governed by the rules 
of quantum mechanics are at best symbolic and at worst misleading. 

where m y = — j, — j + 1, . . . , j — 1, j and where j is the total angular momen- 
tum quantum number. From Equations 2.68, 2.70, and 2.73 it is apparent that 

mj = m,+ m s = m f ± \ (2.74) 

Since is always an integer, must be half-integral ( + ± f , + f , . . . ) and 

thus j must be half-integral. The vector coupling of Equation 2.71 suggests only 
two possible values for j: t + \ or t — which are illustrated in Figure 2.19. 

Usually, we indicate the j value as a subscript in spectroscopic notation. Thus, 
for (= 1 (p states), there are two possible j values: \ = \ and t— \ = We 

would indicate these states as p 3/2 and p 1/2 . When there is an additional 
quantum number, such as a principal quantum number n (or perhaps just an 
index which counts the states in order of increasing energy), we indicate it as 
2p 3/2 , 3p 3/2 , and so on. 

In atoms, it is often useful for us to picture the electrons as moving in well 
defined orbits with definite t and j. It is not at all obvious that a similar picture 
is useful for nucleons inside the nucleus, and thus it is not clear that t and j will 
be useful labels. We discuss this topic in detail when we consider the nuclear shell 
model in Chapter 5. 



2.6 PARITY 

The parity operation causes a reflection of all of the coordinates through the 
origin: r -» -r. In Cartesian coordinates, this means x -> -x, y -» — y , z 
— z; in spherical coordinates, r -> r, 6 -> 7T — 0, <f>-><£ + 77. If a system is left 
unchanged by the parity operation, then we expect that none of the observable 
properties should change as a result of the reflection. Since the values we measure 
for the observable quantities depend on |^| 2 , then we have the following 
reasonable assertion: 

If V(r) = V(~r), then |^(r )| 2 = |i|/(-r) | 2 . 
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This assertion, whose reverse is also true, has two important consequences for our 
work in nuclear physics: 



1. If |t//(r)| 2 = |^( — r)| 2 fh en ^(~ r ) = That is, the parity operation 

has either of two effects on a wave function. The case ^(-r) = + i//(r) is 
known as positive or even parity, while the case \p( — r)= — ip(r) is negative 
or odd parity. If the potential V(r) is left unchanged by the parity operation, 
then the resulting stationary-state wave functions must be of either even or 
odd parity. Mixed-parity wave functions are not permitted. Recall our 
solutions for the one-dimensional harmonic oscillator. The potential \kx 2 is 
certainly invariant with respect to the parity operation x -> — x. The wave 
functions listed in Table 2.1 have either only odd powers of x, and therefore 
odd parity, or only even powers of x, and therefore even parity. Polynomials 
mixing odd and even powers do not occur. Also, review the solutions for the 
finite potential well. Since the well lies between x = +<3/2 and x = —a/2, 
the potential is symmetric with respect to the parity operation: V(x) = 
V( — x). Notice the solutions illustrated in Figure 2.8. For some of the 
solutions, ip( — x) = v//(x) and their parity is even; the other solutions have 
\p(-x) = \p(x) and odd parity. 

In three dimensions, the parity operation applied to the Y € gives a phase 



(-ir 









0 ,/ 



8 ,* + «) = (- 



(2.75) 



Central potentials, which depend only on the magnitude of r, are thus 
invariant with respect to parity, and their wave functions have definite parity, 
odd if £ is odd and even if £ is even. 

The wave function for a system of many particles is formed from the 
product of the wave functions for the individual particles. The parity of the 
combined wave function will be even if the combined wave function repre- 
sents any number of even-parity particles or an even number of odd-parity 
particles; it will be odd if there is an odd number of odd-parity particles. 
Thus nuclear states can be assigned a definite parity, odd or even. This is 
usually indicated along with the total angular momentum for that state, as 
for example, § + or \ In Chapter 10 we will discuss how the parity of a 
state can be determined experimentally. 

2. The second consequence of the parity rule is based on its converse. If we find 
a system for which |^(r)| 2 + |^( — r)| 2 , then we must conclude that V(r ) =£ 
K( — r); that is, the system is not invariant with respect to parity. In 1957 it 
was discovered that certain nuclear processes (fi decays) gave observable 
quantities whose measured values did not respect the parity symmety. On the 
other hand, no evidence has yet been obtained that either the strong nuclear 
interaction or the electromagnetic interaction violate parity. The establish- 
ment of parity violation in ft decay was one of the most dramatic discoveries 
in nuclear physics and has had profound influences on the development of 
theories of fundamental interactions between particles. A description of these 
experiments is given in Section 9.9. 




ELEMENTS OF QUANTUM MECHANICS 39 



2.7 QUANTUM STATISTICS 

When we group several particles together to make a larger quantum system 
(several nucleons in a nucleus, several electrons in an atom, several atoms in a 
molecule) a new quantum effect arises if the particles are indistinguishable from 
one another. Let us consider the case of two particles, for example, the two 
electrons in a helium atom. Suppose one electron is described by coordinates r x 
and is in the state \p A , while the other electron is described by coordinates r 2 and 
is in the state \p B . The combined wave function is the product of the two 
component wave functions; thus ip = 'p A ( r i) i> B (r 2 ). Now suppose the two elec- 
trons are interchanged so that the new wave function is \p' = \p B (r x )ip A (r 2 ). Is 
there any measurement we could do to detect whether this interchange had taken 
place? 

If the electrons are truly indistinguishable , the answer to this question must be 
no. There is no observational scheme for distinguishing the “ first electron” from 
the “second electron.” Thus we have a result that is somewhat similar to our 
result for the parity operation: Probability densities must be invariant with respect 
to exchange of identical particles . That is, the exchanged wave function \p 21 can at 
most differ only in sign from the original wave function xp 12 . We therefore have 
two cases. If the sign does not change upon exchange of the particles, we have a 
symmetric wave function; for symmetric wave functions, \p 12 — \p 2l . If the ex- 
change changes the sign, we have an antisymmetric wave function, for which 
xpii — —*P i 2 * All combined wave functions representing identical particles must be 
either completely symmetric or completely antisymmetric. No “mixed symmetry” 
wave functions are allowed. 

When we turn to our laboratory experiments to verify these assertions, we find 
a further classification to which there are no known exceptions: all particles with 
integral spins (0, 1, 2, ... ) have symmetric combined wave functions, while all 
particles with half-integral spins f , § , . . . ) have antisymmetric combined wave 

functions. 

The above two-particle functions \p and \p' will not do for combined wave 
functions because they are neither symmetric nor antisymmetric. That is, \p' does 
not at all look like either ip or — ip. Instead, consider the following combined 
wave function: 



1 r 

^12 = Ai r l) ^ B^l) i x pA( r 2)] (2-76) 

If we choose the plus sign, then the combined wave function is symmetric with 
respect to interchange of the particles. If we choose the minus sign then the result 
is an antisymmetric wave function. The factor of 1 / fil ensures that the resulting 
combination is normalized (assuming that each of the component wave functions 
is itself normalized). 

A special case arises when we have identical quantum states A and B. (We can 
regard A and B as representing a set of quantum numbers.) When A is the same 
as B , the antisymmetric combination wave function vanishes identically, and so 
its probability density is always zero. The probability to find two identical particles 
of half-integral spin in the same quantum state must always vanish. This is of course 
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just the Pauli exclusion principle , which determines why atomic subshells fill in a 
certain way. This vanishing of the antisymmetric wave function is the mathemati- 
cal basis of the Pauli principle. No such vanishing occurs for the symmetric 
combination, so there is nothing to prevent identical particles of integral spin 
from occupying the same quantum state. 

Later in this text, we apply the Pauli principle to nucleons and show its 
importance in understanding the nuclear shell model. We also construct some 
simple antisymmetric wave functions for the quarks that make up nucleons and 
other similar particles. 



2.8 TRANSITIONS BETWEEN STATES 

A true stationary state lives forever. The expectation values of physical observ- 
ables, computed from the wave function of a stationary state, do not change with 
time. In particular, the expectation value of the energy is constant in time. The 
energy of the state is precisely determined, and the uncertainty in the energy, 

A E = /(£ 2 ) - (E) 2 (2.77) 

vanishes, because ( E 2 ) = ( E ) 2 for this case. The Heisenberg relationship, 
A£A/ > h/ 2, then implies that A / = oo. Thus a state with an exact energy lives 
forever; its lifetime against decay (to lower excited states, for example) is infinite. 

Now suppose our system is subject to a weak perturbing potential V\ in 
addition to the original potential V. In the absence of V\ we can solve the 
Schrodinger equation for the potential V and find a set of eigenstates \p n and 
corresponding eigenvalues E n . If we now include the weak additional potential 
V\ we find that the states are approximately, but not exactly, the previous 
eigenstates \p n of V. This weak additional potential permits the system to make 
transitions between the “approximate” eigenstates xp n . Thus, under the interac- 
tion with a weak electromagnetic field, a hydrogen atom can make transitions, 
such as 2p -> Is or 3d -> 2p. We still describe the various levels as if they were 
eigenstates of the system. 

Even though a system may make a transition from an initial energy state E i to 
a final state E f , energy must be conserved. Thus the total decay energy must be 
constant. If the final state E f is of lower energy than E { , the energy difference 
E { — E { must appear as radiation emitted in the decay. In transitions between 
atomic or nuclear excited states, a photon is emitted to carry the energy E l — E t . 

A nonstationary state has a nonzero energy uncertainty A E. This quantity is 
often called the “width” of the state and is usually represented by T. The lifetime 
t of this state (the mean or average time it lives before making a transition to a 
lower state) can be estimated from the uncertainty principle by associating r with 
the time A t during which we are permitted to carry out a measurement of the 
energy of the state. Thus r — h/Y. The decay probability or transition probability 
A (the number of decays per uirirtime) is inversely related to the mean lifetime t: 

1 

A = - (2.78) 

T 

It would be very useful to have a way to calculate A or t directly from the 
nuclear wave functions. We can do this if we have knowledge of (1) the initial 
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and final wave functions and \p f , which we regard as approximate stationary 
states of the potential V; and (2) the interaction V' that causes the transition 
between the states. The calculation of A is too detailed for this text, but can be 
found in any advanced text on quantum mechanics. We will merely state the 
result, which is known as Fermi's Golden Rule : 

A= T |F fi '| 2 p(£ f ) (2.79) 

The quantity V { - has the form of an expectation value: 

= / 'P* V'xf/j dv (2.80) 

Notice again the ordering of the states f and i in the integral. The integral V { [ is 
sometimes called the matrix element of the transition operator V'. This terminol- 
ogy comes from an alternative formulation of quantum mechanics based on 
matrices instead of differential equations. Be sure to take special notice that the 
decay probability depends on the square of the transition matrix element. 

The quantity p(E ( ) is known as the density of final states. It is the number of 
states per unit energy interval at E ( , and it must be included for the following 
reason: if the final state E ( is a single isolated state, then the decay probability 
will be much smaller than it would be in the case that there are many, many 
states in a narrow band near E t . If there is a large density of states near £ f , there 
are more possible final states that can be reached by the transition and thus a 
larger transition probability. The density of final states must be computed based 
on the type of decay that occurs, and we shall consider examples when we discuss 
/? decay, y decay, and scattering cross sections. 
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PROBLEMS 

1. Derive Equation 2.37 and plot the transmission coefficient as a function of 
the energy E of the incident particle. Comment on the behavior of T. 

2 . Derive Equation 2.39 and plot the transmission coefficient as a function 
of E. 

3. Solve the Schrodinger equation for the following potential: 

V(x) =oo x <0 

= — V 0 0 <x < a 
= 0 x > a 

Here V 0 is positive and solutions are needed for energies E > 0. Evaluate all 
undetermined coefficients in terms of a single common coefficient, but do 
not attempt to normalize the wave function. Assume particles are incident 
from x — — oo . 

4 . Find the number of bound states and their energies in the finite one-dimen- 
sional square well when P = 10. 

5. Find the solution to the “half” harmonic oscillator: 

V(x) = oo x < 0 

= \kx 2 x > 0 

Compare the energy values and wave functions with those of the full 
harmonic oscillator. Why are some of the full solutions present and some 
missing in the “half” problem? 

6. For the ground state and first two excited states of the one-dimensional 
simple harmonic oscillator, find the probability for the particle to be beyond 
the classical turning points. 

7 . (a) For the ground state of the one-dimensional simple harmonic oscillator, 

evaluate (x) and (x 2 ). 

(b) Find Ax = [<x 2 > - <x> 2 ] 1/2 . 

(c) Without carrying out any additional calculations, evaluate (pS) and 
</>*>• ( Hint: Find {p 2 J2m).) 

(d) Evaluate A p x and the product Ax • A/? x . A wave packet with this shape 
(called a Gaussian shape) is known as a “minimum-uncertainty” wave 
packet. Why? 

8. (a) Find the wave functions and energy levels for a particle confined to a 

two-dimensional rectangular box, with 

V(x, y) = 0 — a < x < + a, —b < y < +b 

= 00 | x | > a , \y I > b 

(b) Make a diagram similar to Figure 2.10 showing the levels and degener- 
acies for b = a and for b = 2a. 

9 . Continue Figure 2.10 upward to 50£ 0 . 

10 . Continue Figure 2.11 upward to 20 E 0 . 
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11 . Carry out the separation procedure necessary to obtain the solution to 
Equation 2.54. 

12 . Show that the first four radial wave functions listed in Table 2.4 are 
solutions to the Schrodinger equation corresponding to the proper value of 
the energy, and also show that they are normalized. 

13 . Find the solutions to the one-dimensional infinite square well when the 
potential extends from —a/2 to +a/2 instead of from 0 to +a. Is the 
potential invariant with respect to parity? Are the wave functions? Discuss 
the assignment of odd and even parity to the solutions. 

14 . Find the angle between the angular momentum vector t and the 2 axis for 
all possible orientations when *?= 3. 

15 . (a) What are the possible values of j for f states? 

(b) What are the corresponding m .? 

(c) How many total states are there? 

(d) How many states would there be if we instead used the labels and 

mp. 

16 . A combined wave function representing the total spin S of three electrons 
can be represented as \p = ^i(f n s )^ 2 ( m s) x p 3 ( m s ) where m s = ± \ for a 
spin- \ electron, (a) List all such possible wave functions and their total 
projection M s . (b) Identify the number of M s values with the number of 
different values of the total spin S. (Some values of S may occur more than 
once.) (c) Draw simple vector diagrams showing how the different couplings 
of s v s 2 , and s 3 can lead to the deduced number and values of S . (Hint: 
First couple two of the spins, and then couple the third to the possible 
resultants of the first two.) (d) In the coupling of four electrons, give the 
possible values of S and their multiplicity, and show that the number of M s 
states agrees with what would be expected from tabulating the possible 
wave functions. 
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NUCLEAR PROPERTIES 



Like many systems governed by the laws of quantum mechanics, the nucleus is a 
somewhat enigmatic and mysterious object whose properties are much more 
difficult to characterize than are those of macroscopic objects. The list of 
instructions needed to build an exact replica of a French colonial house or a ’57 
Chevy is relatively short; the list of instructions necessary to characterize all of 
the mutual interactions of the 50 nucleons in a medium weight nucleus could 
contain as many as 50! or about 10 64 terms! We must therefore select a different 
approach and try to specify the overall characteristics of the entire nucleus. Are 
there a few physical properties that can be listed to give an adequate description 
of any nucleus? 

To a considerable extent, we can describe a nucleus by a relatively small 
number of parameters: electric charge, radius, mass, binding energy, angular 
momentum, parity, magnetic dipole and electric quadrupole moments, and 
energies of excited states. These are the static properties of nuclei that we 
consider in this chapter. In later chapters we discuss the dynamic properties of 
nuclei, including decay and reaction probabilities. To understand these static and 
dynamic properties in terms of the interaction between individual nucleons is the 
formidable challenge that confronts the nuclear physicist. 



3.1 THE NUCLEAR RADIUS 

Like the radius of an atom, the radius of a nucleus is not a precisely defined 
quantity; neither atoms nor nuclei are solid spheres with abrupt boundaries. Both 
the Coulomb potential that binds the atom and the resulting electronic charge 
distribution extend to infinity, although both become negligibly small at distances 
far beyond the atomic radius (10 ~ 10 m). What is required is an “operational 
definition” of what we are to take as the value of the atomic radius. For instance, 
we might define the atomic radius to be the largest mean radius of the various 
electronic states populated in an atom. Such a property would be exceedingly 
difficult to measure, and so more practical definitions are used, such as the 
spacing between atoms in ionic compounds containing the atom of interest. This 
also leads to some difficulties since we obtain different radii for an atom when it 
is in different compounds or in different valence states. 
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For nuclei, the situation is better in some aspects and worse in others. As we 
will soon discuss, the density of nucleons and the nuclear potential have a similar 
spatial dependence — relatively constant over short distances beyond which they 
drop rapidly to zero. It is therefore relatively natural to characterize the nuclear 
shape with two parameters: the mean radius, where the density is half its central 
value, and the “skin thickness,” over which the density drops from near its 
maximum to near its minimum. (In Section 5 of this chapter, we discuss a third 
parameter, which is necessary to characterize nuclei whose shape is not spherical.) 

The problems we face result from the difficulty in determining just what it is 
that the distribution is describing; the radius that we measure depends on the 
kind of experiment we are doing to measure the nuclear shape. In some 
experiments, such as high-energy electron scattering, muonic X rays, optical and 
X-ray isotope shifts, and energy differences of mirror nuclei, we measure the 
Coulomb interaction of a charged particle with the nucleus. These experiments 
would then determine the distribution of nuclear charge (primarily the distribution 
of protons but also involving somewhat the distribution of neutrons, because of 
their internal structure). In other experiments, such as Rutherford scattering, a 
decay, and pionic X rays, we measure the strong nuclear interaction of nuclear 
particles, and we would determine the distribution of nucleons, called the 
distribution of nuclear matter . 

The Distribution of Nuclear Charge 

Our usual means for determining the size and shape of an object is to examine 
the radiation scattered from it (which is, after all, what we do when we look at an 
object or take its photograph). To see the object and its details, the wavelength of 
the radiation must be smaller than the dimensions of the object; otherwise the 
effects of diffraction will partially or completely obscure the image. For nuclei, 
with a diameter of about 10 fm, we require A < 10 fm, corresponding to p > 100 
MeV/c. Beams of electrons with energies 100 MeV to 1 GeV can be produced 
with high-energy accelerators, such as the Stanford linear accelerator, and can be 
analyzed with a precise spectrometer to select only those electrons that are 
elastically scattered from the selected nuclear target. Figure 3.1 shows an example 
of the results of such an experiment. The first minimum in the diffractionlike 
pattern can clearly be seen; for diffraction by a circular disk of diameter D , the 
first minimum should appear at 6 = sin -1 (1.22 A /D), and the resulting estimates 
for the nuclear radii are 2.6 fm for 16 0 and 2.3 fm for 12 C. These are, however, 
only rough estimates because potential scattering is a three-dimensional problem 
only approximately related to diffraction by a two-dimensional disk. 

Figure 3.2 shows the results of elastic scattering from a heavy nucleus, 208 Pb. 
Several minima in the diffractionlike pattern can be seen. These minima do not 
fall to zero like diffraction minima seen with light incident on an opaque disk, 
because the nucleus does not have a sharp boundary. 

Let us try to make this problem more quantitative. The initial electron wave 
function is of the form e lkx ' r , appropriate for a free particle of momentum 
p x = hk v The scattered electron can also be regarded as a free particle of 
momentum = hk t and wave function e lkrr . The interaction V(r) converts the 
initial wave into the scattered wave, and according to Equation 2.80 the probabil- 
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Figure 3.1 Electron scattering from ie O and 12 C. The shape of the cross section 
is somewhat similar to that of diffraction patterns obtained with light waves. The 
data come from early experiments at the Stanford Linear Accelerator Center (H. F. 
Ehrenberg et al., Phys. Rev. 113 , 666 (1959)). 

ity for the transition will be proportional to the square of the quantity 

F{k { , k f ) = Jipf V(r )ip i dv (3.1) 



F(q) = j e“ l ' r V(r) dv 



(3.2) 



apart from a normalization constant, which is chosen so that F(0) = 1. Here 
q = k x — /r f , which is essentially the momentum change of the scattered electron. 
The interaction V{r) depends on the nuclear charge density Zep e (r '), where r' is 
a coordinate describing a point in the nuclear volume and p e gives the distribu- 
tion of nuclear charge. That is, as indicated in Figure 3.3, an electron located at r 
feels a potential energy due to the element of charge dQ located at r': 

edQ 

dV — 

4w ok - ''I 

Ze 2 p e (r') dv' 

4ire 0 \r - r'\ 

To find the complete interaction energy V(r), we sum over all of the contribu- 
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Figure 3.2 Elastic scattering of electrons from 208 Pb. Note the different vertical 
and horizontal scales for the two energies. This also shows diffractionlike behavior, 
but lacks sharp minima. (J. Heisenberg et al., Phys. Rev. Lett. 23, 1402 (1969).) 



tions dQ within the nucleus: 

V{r) 



Ze 2 r Pe(r')dv' 

4 ve 0 J |r-r'| 



(3.4) 



Writing q*r= qr sin 8 in Equation 3.2 and integrating over r, the properly 
normalized result is 



F(q) = f e*- r 'p e {r’)dv’ 



(3.5) 



and if p e (r ') depends only on the magnitude r' (not on 6' or <j>') we obtain 

4?r r / 

F(q ) = — sin qr'p e (r')r' dr' (3.6) 

q J 

This quantity is a function of q , the magnitude of q. Since we originally assumed 
that the scattering was elastic, then \p { \ = |p f | and q is merely a function of the 
scattering angle a between p { and p f ; a bit of vector manipulation shows 

q = ( 2 p/h) sin a/2 where p is the momentum of the electron. Measuring the 

scattering probability as a function of the angle a then gives us the dependence 
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Figure 3.3 The geometry of scattering experiments. The origin of coordinates is 
located arbitrarily. The vector r' locates an element of charge dQ within the 
nucleus, and the vector r defines the position of the electron. 

of Equation 3.6 on q. The quantity F(q) is known as a form factor, and the 
numerical inversion of Equation 3.6, actually an inverse Fourier transformation, 
then gives us p e (r'). 

The results of this procedure for several different nuclei are shown in Figure 
3.4. One remarkable conclusion is obvious — the central nuclear charge density is 
nearly the same for all nuclei. Nucleons do not seem to congregate near the 
center of the nucleus, but instead have a fairly constant distribution out to the 
surface. (The conclusion from measurements of the nuclear matter distribution is 
the same.) Thus the number of nucleons per unit volume is roughly constant : 

A 

y - constant (3.7) 

j77 R 

where R is the mean nuclear radius. Thus R cc A 1/3 , and defining the proportion- 
ality constant R 0 gives 

R = R 0 A 1/3 (3.8) 

From electron scattering measurements, such as those in Figure 3.4, it is 
concluded that R 0 - 1.2 fm. These measurements give the most detailed descrip- 
tions of the complete nuclear charge distribution. 

Figure 3.4 also shows how diffuse the nuclear surface appears to be. The charge 
density is roughly constant out to a certain point and then drops relatively slowly 
to zero. The distance over which this drop occurs is nearly independent of the 
size of the nucleus, and is usually taken to be constant. We define the skin 
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Figure 3.4 The radial charge distribution of several nuclei determined from 
electron scattering. The skin thickness t is shown for O, Ni, and Pb; its value is 
roughly constant at 2.3 fm. The central density changes very little from the lightest 
nuclei to the heaviest. These distributions were adapted from R. C. Barrett and D. F. 
Jackson, Nuclear Sizes and Structure (Oxford: Clarendon, 1977), which gives more 
detail on methods of determining p(r). 



thickness parameter t as the distance over which the charge density falls from 
90% of its central value to 10%. The value of t is approximately 2.3 fm. 

Figure 3.5 shows a more quantitative determination of the relationship be- 
tween the nuclear radius and mass number, based on electron scattering results. 
The root mean square (rms) radius, (r 2 ) 1/2 , is deduced directly from the 
distribution of scattered electrons; for a uniformly charged sphere ( r 2 ) = f R 2 , 
where R is the radius of the sphere. Figure 3.5 shows that the dependence of R 
on A l/3 is approximately valid over the range from the lightest to the heaviest 
nuclei. From the slope of the line, we deduce R 0 = 1.23 fm. 

The nuclear charge density can also be examined by a careful study of atomic 
transitions. In solving the Schrodinger equation for the case of the atom with a 
single electron, it is always assumed that the electron feels the Coulomb attrac- 
tion of a point nucleus, V(r) = — Ze 2 /47re 0 r. Since real nuclei are not points, 
the electron wave function can penetrate to r < R, and thus the electron spends 
part of its time inside the nuclear charge distribution, where it feels a very 
different interaction. In particular, as r 0, V(r) will not tend toward infinity 
for a nucleus with a nonzero radius. As a rough approximation, we can assume 
the nucleus to be a uniformly charged sphere of radius R , for which the potential 
energy of the electron for r < R is 



V'{r) 



Ze 2 f 3 
4vt£ 0 « \ 2 



2 \R 



(3.9) 



while for r > R, the potential energy has the point-nucleus form. 

The total energy E of the electron in a state of a point nucleus depends in 
part on the expectation value of the potential energy 



(V) = j\p* V\p n dv 



(3.10) 



